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Abstract 

The residual mass of the pion in a finite spatial box at vanishing quark masses, the mass gap, 
is computed with two flavors of dynamical clover fermions. The result is compared with 
predictions of chiral perturbation theory in the 8 regime. 



Understanding finite size effects has always been an important issue in lattice studies of QCD. 
For an accurate determination of hadron observables it is of foremost importance to account for 
finite volume corrections. In addition, finite size effects can provide valuable information on the 
low-energy physics of the system. As lattice calculations are reaching the physical pion mass, 
the impact of finite volume will become pronounced. The main effect is caused by modified pion 
dynamics arising from the boundary conditions being imposed on the system. 
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In a finite spatial box chiral symmetry does not break down spontaneously. This results in 
a mass gap, even if the quark masses vanish. The physics behind that is described by a simple 
quantum mechanical rotator JT], whose energy levels can be computed from an expansion of the 
chiral effective theory in the 6 regime m n L <sc 1 , T » L, where L (T) is the spatial (temporal) extent 
of the box. We are now in the position to probe the pion mass near the chiral limit and account for 
this effect. What makes this problem attractive, beyond the computational task of solving QCD in 
a finite volume, is that it is a universal feature of quantum mechanics, which is only constrained 
by the symmetry of the system. For a previous attempt of extracting the mass gap see 0. 

In a series of papers Leutwyler (T), Hasenfratz and Niedermayer []3l and Hasenfratz [4] have 
computed the energy levels for two flavors of dynamical quarks to leading (L), next-to-leading 
(NL) and next-to-next-to-leading (NNL) order. The mass gap, i.e. the residual mass of the pion in 
the chiral limit, up to NNL order turns out to be 
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where F n is the pion decay constant, and A, are the intrinsic scale parameters of the low-energy 
constants 0, 

Our simulations are performed with the Wilson gauge action and two flavors of mass-degenerate 
0(a) improved Wilson fermions on a range of lattice volumes with periodic (antiperiodic) bound- 
ary conditions in the spatial (temporal) direction. Two values of the gauge coupling, J3 = 5.29 and 
5.40, have been selected here. The hopping parameters and the lattice volumes of the individual 
runs are listed in Tabled] It is convenient to express the scale in terms of the force parameter r . 
In the chiral limit we find r /a = 6.201(25) and 6.946(44), respectively. We set the scale by the 
nucleon mass, which gives r = 0.467 fm [6]|. This translates into lattice spacings a = 0.075 and 
0.067 fm. 

The quark masses are computed from the ratio (see e.g. 0, and references therein) 
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at t » a, where ^ is the improved axial vector current (HI, 
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Table 1: Parameters of the lattice data sets used in this analysis. 



and and P are the standard axial vector current and pseudoscalar density, respectively. The 
quark masses show practically no finite size effects beyond the statistical errors. In Table[2]we give 
the results from the largest lattices. Unlike the quark masses, the pion masses show significant 
finite size effects. The results obtained on the various volumes are shown in Table |2] as well. The 
pion masses from the lattices in Table[T|are printed in roman. They reach as low as 170 MeV. 

Some of the pion masses at the larger quark masses are only known on 16 3 and 24 3 lattices. 
Furthermore, we miss one pion mass at the smallest quark mass on the 24 3 lattice at /3 = 5.40. 
In these cases the pion masses have been extrapolated to the 'missing' volumes by means of the 
0(p 4 ) finite size shift formula [|9l 



valid in the p regime m n L » 1 , with 
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Table 2: The pion and bare quark masses on the various lattices. The roman numbers are from our 
simulations on the quoted volumes, while the italic numbers have been obtained by extrapolation 
as described in the text. 
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5°(i) = 2K X (A) , 5 2 (i) = 2K 2 (A)jA . (9) 
We take F n = 92.4 MeV throughout this paper. The Taylor coefficients g, are 

n n 1 1 7r 

8o = 2--, gi = - - - , ^2 = 2~g- ( 10 ) 

The first (second) expression in the square brackets of eq. © corresponds to 0(p 2 ) (0(p 4 )) cor- 
rections in the Lagrangian. The low-energy parameters l\, I 2 and / 4 are taken from Q: 

h^-0A, Z 2 ^4.3, l 4 ^4A, (11) 
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Figure 1: The pion mass squared against the quark mass on the 24 3 lattice at /S = 5.40, together 
with the chiral fit (fT2j) . The circles refer to the roman numbers in Table [2l obtained directly on the 
quoted lattice, while the square is the extrapolated value and refers to the italic number in Tabled 




Figure 2: The pion mass squared against the quark mass on the 24 3 lattice at /S = 5.29, together 
with the chiral fit (fT2~l) . The circles refer to the roman numbers in Table [2l obtained directly on the 
quoted lattice. The extrapolated value referring to the italic number in Table [2] is not shown here. 



5 




0.1 0.2 0.3 0.4 



r m q 

Figure 3: The pion mass squared against the quark mass on the 32 3 lattice at [S = 5.40, together 
with the chiral fit (fT2j) . The circles refer to the roman numbers in Table El obtained directly on 
the quoted lattice, while the squares are the extrapolated values and refer to the italic numbers in 
Tabled 
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Figure 4: The pion mass squared against the quark mass on the 40 3 lattice at /S = 5.29, together 
with the chiral fit (fT2j) . The circles refer to the roman numbers in Tabled obtained directly on 
the quoted lattice, while the squares are the extrapolated values and refer to the italic numbers in 
Tabled 
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Figure 5: The mass gap as a function of lattice size for (8, L/a)=(5A0,24), (5.29,24) (5.40,32) 
(5.29,40), from left to right. The solid (dashed) [dotted] curve is the prediction of the chiral 
effective theory to NNL (NL) [L] order. The shaded area corresponds to m n L < 1, indicating the 6 
regime. The dashed-dotted horizontal line marks the position of the physical pion mass. We have 
set F n = 92.4 MeV, and the scale parameters Ai, A 2 have been taken from [5], as stated in eq. (fTTI) . 

while Z 3 = 4.2 is taken from our fit of the pion masses in Figs. |T|— [4]( see Table©. The extrapolated 
numbers are printed in italics in Table |2] As the extrapolation concerns mainly the larger quark 
masses, any small error inherent in the extrapolation should not affect our final result significantly. 
All masses in question are in the p regime. 

In Figs. [T|-|4] we pl ot me pi° n ma ss squared against the quark mass on the 24 3 lattices at 
8 = 5.40 and 8 = 5.29, on the 32 3 lattice at 8 = 5.40, and on the 40 3 lattice at 8 = 5.29. We fit the 
data by the formula 

(r m„) 2 = A + B r m q [l + C r Q m q ln(D r m q )\ , (12) 

which allows for a residual pion mass, A = (r m™ s ) 2 , in the chiral limit, and at larger quark masses 
turns over to the 0(p A ) chiral expansion valid in the p regime Q. This is to say, eq. (fT2l) interpolates 
between the energy levels in the 8 and the p regime. 

As we are working at fixed L, we are actually fitting the mass dependence of m n (L)\ L rather 
than m n , and one might expect corrections to the fit formula, which cannot be accounted for by our 
ansatz (fT2l) . However, this is not the case. In Fig.[6]in the Appendix we compare the leading finite 
size correction in the p regime with the leading contribution to the mass gap in the 6 regime. At 
our smallest quark mass the correction to m n (in fact to the difference m n - m™ s ) turns out to be a 
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Table 3: The residual pion mass m^ es and the low-energy constant / 3 . 

fraction of the statistical error only for L/a = 40 and 32, and about the size of the statistical error 
for L/a = 24. Moreover, the corrections in the p regime are about a factor of 12 smaller than the 
mass gap. 

The individual fits are shown in Figs. [Q-IU On each of our lattices we obtain a residual mass, 
which is distinctly above zero. The results are summarized in Table [31 Barring any finite size 
effects, the low-energy constant Z 3 is given by 

/_ 3 = ln (4^) ' ln(roA 3 ) 2 = 327r 2 (r OJ F ;r ) 2 ^ln(|). (13) 

The numbers of the individual fits are listed in Table [31 The average value turns out to be h = 
4.2(2), which lies at the upper end of the results reported in the literature IfTOl . 

In Fig. [5] we finally compare the residual pion masses of Table [3] with the prediction (OQ) of the 
chiral effective theory. We notice that the chiral series seems to have converged for L/r > 5, and 
below that the NL and NNL order results differ only by a few percent. We find good agreement 
with the NNL order formula flU for all our lattice volumes, ranging from L/r « 3.5 to « 6.5. 
This is a bit surprising, as our results are barely touching the 6 regime. The horizontal line in 
Fig. \5\ shows the position of the physical pion mass. To capture the physics of light pions, and the 
spontaneous breakdown of chiral symmetry, the spatial extent of the lattice should be L > 3 fm, 
such that «: m£ hys . 

With precise data at (say) L/r w 6, and quark masses at - or close to - the physical point, it 
should be possible to determine F n accurately. The advantage of such a calculation is that it does 
not demand renormalization of the axial vector current. A fit of © to m^ es on our largest lattice, 
L/a = 40 at = 5.29 (the rightmost data point in Fig. [5]), gives the pion decay constant in the 
chiral limit 

F = F„\ miT=0 = 78±Jo MeV (14) 

(still using r = 0.467 fm), which is in the right ball-park IfTTII . but not accurate enough to be useful 
yet. Due to the relatively small correction of the NNL order contribution it will, however, not be 
possible to determine the low-energy constants l\ and l 2 reliably. 



8 




Figure 6: The leading contribution to the finite size mass shift of the pion, in units of (F^L 3 ) -1 , as 
a function of m n L (solid and dashed curve), compared with the leading order contribution to the 
mass gap (solid circle). Within the shaded area m„(L) x F 2 l} drops from 3/2 at m n L = to 1/8 at 
\IF 2 n L 2 « m^L « 1. 
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Appendix 



In Fig. [6] we show the leading 0(p 2 ) finite size correction to the mass of the pion, 



n fe m i 
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together with the mass gap (or energy gap) to leading order, 

< ! ~. (16) 

According to eq. (25) of [Q]|, the finite volume shift drops by a factor of 12 from 3/2F^L 3 to 
1 /8F^L 3 in the shaded region. Incidentally (?), the p regime finite size correction formula (fT5l) 
extrapolates to exactly the same value, 1/8F^L 3 , in the chiral limit. For small m n L we can express 
the sum over n in eq. (fT5l ) by an integral: 

m n {L) = m n + dvAnv - — . (17) 

\6n z F^ Jo m n vh 

Changing the variable v to p = m n vL, we obtain 

1 f°° 

ra^L) = m n + 4nF 2 L 3 J d/j.fiKi(pi). (18) 
The integral is known analytically, 

dfi/uK l ( M ) = ^, (19) 

giving 

m n {L) = m n + —— (20) 
8F2L 3 

as the small m n L limit of (fT5l) . The 0(p 4 ) contribution to the mass shift would be suppressed by 
another factor of L' 2 . 

An analytical curve, that smoothly connects the mass gap of the 6 regime with the 0(p 4 ) cor- 
rection formula of the p regime, would be very valuable. 
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